Pis'ma v ZhETF 



A search for Extragalactic Sources of Ultrahigh-Energy Cosmic Rays 
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Possible extragalactic sources of cosmic rays at energies above 4 x 10 eV detected with the Yakutsk array 
are sought. Correlation of the shower arrival directions with objects from Veron's catalog that are located 
closer than 100 Mpc from the Earth confirms the observations at the Pierre Auger observatory, as well as the 
Greisen-Zatsepin-Kuzmin effect on ultrahigh-energy cosmic rays. The detailed analysis of the data reveals the 
classes of objects belonging to the active galactic nuclei that are probable sources of ultrahigh-energy cosmic 
rays. 

PACS: 96.40.Pq, 98.54.Cm, 98.70.Sa 



The problem of the origin of ultrahigh-energy cos- 
mic rays (UHECRs) remains unsolved in spite of con- 
siderable efforts of researchers. The reason is that the 
observed distribution of the arrival directions of cos- 
mic rays is practically isotropic, due to the trajectories 
of charged particles deflected in magnetic helds in the 
Galaxy and beyond. Indications of a possible corre- 
lation of UHECR arrival directions with certain most 
active extragalactic objects have been obtained only at 
the highest energies E > 4 x 10 19 eV (= 40 EeV), where 
the deflection of protons, the most probable particles of 
cosmic rays, is comparable or smaller than the angular 
resolution of detectors. 

In particular, Farrar and Biermann [1] found corre- 
lation of the highest energy (E > 10 199 eV) extensive 
,air showers (EASs) detected to that time with com- 
pact radio quasars at high rcdshifts. On the contrary, 
Glushkov [2] pointed out the correlations with closer 
quasars, z < 2.5. However, these results were not con- 
firmed by subsequent observations. 

Another type of the extragalactic objects that can 
generate UHECRs is presented by gamma ray bursts [3] . 
In this search for correlation is complicated, 

because these are short-term events with duration less 
than several minutes. Nevertheless, the data from the 
Pierre Auger observatory (PAO, 609 161 EASs) were 
compared 4J with the observations by Swift, HETE, 
INTEGRAL, etc. missions (284 bursts detected with 
an angular resolution better than 1°). To reveal corre- 
lation of the arrival directions with the coordinates of 
the bursts within the angular range from 5° to 30°, a 
time interval of 100 days before and after each burst 
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is taken. Comparing the event rate of cosmic rays be- 
fore and after the birst as a function of time difference, 
Anchordoqui [3] concluded that there is no frequency 
difference detectable. 

Stanev et al. [5] noted that the existing set of cos- 
mic rays with energies E > 40 EeV detected before 1995 
tends to align along the Supergalactic plane, where the 
density of galaxies is relatively high. However, analy- 
sis of data from five Northern-Hemisphere arrays (114 
showers, E > 40 EeV) performed in 2000 [B] revealed 
no significant excess of particles from the plane. At the 
same time, Uchihori et al. [BJ noted that the number 
of doublets and triplets (coincidences in the arrival di- 
rections of two and three particles, respectively, within 
4°) in ±10° band near the Supergalactic plane is larger 
than that expected for the isotropic distribution. They 
treated it as the possibility for the part of cosmic rays 
at energies above 4 x 10 19 eV to correlate with the Su- 
pergalactic plane. 

Blasars, including BL Lacertae (BL Lacs) and OVV 
quasars, used to be regarded as active galactic nuclei 
(AGNs) with relativistic jets pointing at the Earth. 
Tinyakov and Tkachev [7 found a significant angular 
correlation of UHECRs detected by the AGASA (39 
EASs, E > 48 EeV) and Yakutsk (26 EASs, E > 24 
EeV) arrays with BL Lacs (m < 18) from catalog [5]. 
The analysis of HiRes data performed by Abbasi et 
al. [5] did not confirm this result for the BL objects, but 
they found a significant correlation (revealed by Gor- 
bunov et al. [10]) of cosmic rays with energies above 
10 19 eV and all events without selection in energy with 
high-polarization (HP) objects, as well as with the com- 
bined BL+HP sample from the same catalog. Abbasi et 
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FigtU Distribution of UHECR arrival directions (cir- 
cles) and active nuclei (squares) in the Galactic coordi- 
nates. The dashed line is the Supergalactic plane and 
the dash-dotted line is bordering the observation region 
of the Yakutsk array. 



al. [5] interpreted this correlation as an indication of the 
possible flux of neutral particles generated by Lacertae. 

Other AGNs, namely Seifert galaxies closer than 40 
Mpc, were treated as possible sources of UHECRs by 
Uryson [IT] . She found angular correlation of the show- 
ers of energies 3.2 x 10 19 < E < 3 x 10 20 eV detected by 
the Akeno, AGASA, Haverah Park, and Yakutsk arrays 
with Seifert galaxies at redshifts z < 0.0092, which are 
faint x-ray and radio sources. Subsequently, correlation 
of UHECR arrival directions with these objects, as well 
as with other possible sources, was re-analyzed in [12] . 
It was shown that correlation with Seifert galaxies is 
observed only in the AGASA data, and if the deflection 
of particles in the Galactic magnetic field is taken into 
account, i.e., under the conditions different from those 
accepted in [TT] . 

The Pierre Auger collaboration (PAC) [13] recently 
analyzed a sample consisting of 81 EASs with energies 
above 40 EeV detected from January 1, 2004 to Au- 
gust 31, 2007. The authors used a part of the data (to 
May 27, 2006) in order to determine the parameters re- 
sulting in the maximum correlation of UHECR arrival 
directions with AGNs. Then, the second part of the 
data was used to confirm the hypothesis obtained. 

As a result, the observed UHECR arrival direc- 
tions are found to be anisotropic, and there is a signif- 
icant correlation of EASs with energies above 56 EeV 
within an angle of ip = 3.1° with AGNs from cata- 
log [8] located at distances z < 0.018 from the Earth 
(closer than 75 Mpc if the Hubble constant is equal to 
71 kms~ 1 Mpc~ 1 ). In the second part of the data (from 
May 27, 2006), 8 of 13 EASs correlate with AGNs under 
the same conditions that have been found for the first 
part of the data, while the number of expected coinci- 
dences is 2.7 in the isotropic case. This corresponds to 
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Fig[2] Excess of UHECRs from the Supergalactic plane 
as a function of energy. The vertical bars are statisti- 
cal errors and the horizontal bars are energy bins. The 
dashed line shows the value expected for the isotropic 
distribution. 

the chance probability P = 1.7 x 10~ 3 for the uniform 
distribution. 

In this work, in order to confirm or reject the PAC 
result assumed as the null hypothesis, an independent 
sample of 51 EASs with energies above 40 EeV and 
zenith angles below 60° that were detected with the 
Yakutsk array and reported by Pravdin et al. Q3] is 
analyzed. The error in the determination of the arrival 
angles of these showers is less than 5°. The energy of a 
primary particle initiating EAS is estimated on the basis 
of the total flux measurement of Cherenkov light and the 
numbers of electrons and muons at the observation level. 
The procedure used was described in [15] [16] [17] . The 
error in the estimation of the energy, SE/E, is about 
30% and less than 50% for the showers with axes in- 
side the array area and in the effective region outside, 
respectively. 

The Galactic coordinates of AGNs from [8] at the 
distances z < 0.015 are shown in Fig. Q]in the Hammer- 
Aitoff projection along with the arrival angles of 22 
EASs with energies above 60 EeV. The observation re- 
gions of the Pierre Auger observatory and Yakutsk ar- 
ray are directed towards and outwards the center of 
the Galaxy, respectively, because these two arrays are 
located in different hemispheres. Thus, these two ob- 
servation regions are complementary and, in particular, 
cover different segments of the Supergalactic plane. The 
showers detected with the Yakutsk array seem to gather 
around the Supergalactic plane: the density of UHECRs 
is higher near this plane. 

To verify this hypothesis by Stanev et al. [5] us- 
ing the extended Yakutsk array database, the excess 
R = (nscp - n ot her)/(nsGP + n ther) of UHECRs ar- 
riving from the ±10° vicinity of the Supergalactic plane, 
tt-sgp, over the number from all other directions, 

Toother i 

is calculated and compared with the value Rq — —0.555 
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FigGD Ratio of the coincidences of UHECR arrival di- 
rections with extragalactic objects to the number of ran- 
dom coincidences expected for the isotropic distribution 
versus the threshold energy of the particles, E > Eth- 
The statistical error bars are also shown. 
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FiglU Ratio of the coincidences of UHECR arrival di- 
rections with AGNs to the expected number in various 
redshift bins, z. The vertical bars are statistical errors. 
The boundaries of z (shown by the horizontal bars) are 
0.001, 0.015, 0.025, 0.05, 0.075, 0.1, 0.15, 0.2, 0.3, 0.4, 
0.5, 0.7, 0.9, 1.5, and 5.4. 



obtained for the isotropic distribution of cosmic rays 
taking into account different exposures of celestial re- 
gions over the array (see Fig. [2]). An exposure was 
calculated using the algorithm described in [IS]. As it 
is seen in the figure, R tends to increase with energy, 
but there is no statistically significant excess of cosmic 
ray flux from the Supergalactic plane in the data in 
any energy bin. Glushkov and Pravdin [19j previously 
claimed that the observed UHECR flux from this plane 
noticeably exceeds the flux expected for the isotropic 
distribution if the angular interval width near the plane 
and energy threshold for the sample are adjusted. In 
particular, the excess in a narrow angular bin (1° — 2°) 
and for the particle energies above 8 EeV is (4 — 5) a for 
the Yakutsk array data. However, the correction fac- 
tor to this excess that should be introduced in order to 
take into account a posteriori optimization of the data 
selection criteria is indefinite in this case. 

Let us proceed to the correlation of UHECRs de- 
tected in Yakutsk P3] with extragalactic objects. The 
test of the PAC hypothesis with the recommended pa- 
rameters (E > 56 EeV, z < 0.018, and ip = 3.1°) shows 
that the arrival directions of 12 of 24 EASs correlate 
with the AGNs from [SJ, while the number of coinci- 
dences expected for the isotropic distribution is 5.6. In 
this case, the chance probability of 12 or more coinci- 
dences is P = 4 x 10~ 3 . Therefore, the Yakutsk array 
data confirm the result obtained by the Pierre Auger 
collaboration, but at a lower significance level. 

Due to the different observation region, as well as the 
energy/ arrival angles estimation procedures, the 'opti- 
mal' parameters of correlation can be different for the 
Yakutsk array data. For this reason, scanning in the 
energy (E > 40 EeV), redshift (0.001 < z < 0.03), 



and angular distance (1° < ip < 6°) is performed to 
determine the maximum ratio of the difference in the 
observed number of coincidences and the number ex- 
pected for the isotropic case to the standard deviation. 
The maximum ratio appears to be reached for 22 EASs 
with energies above 60 EeV, 12 of which arrive within 
-0 = 3° of the AGNs (while an expected number is 4.1) 
at distance from the Earth less than z — 0.015 (63 Mpc). 
The chance probability is P = 2 x 10~ 4 . In this case, 
it is also necessary to determine a penalty factor to the 
probability due to a posteriori selection of the parame- 
ters; or the chosen parameters should be used as a null 
hypothesis for the testing with independent data. 

Figure [3] shows the ratio of the observed number 
of coincidences of UHECR arrival directions (hereafter, 
within 3°) with quasars, Lacertae, and AGNs from [SJ 
to the number of random coincidences expected for the 
isotropic distribution, taking into account the exposure 
of sky regions to the Yakutsk array. The active nuclei 
were chosen at the distances z < 0.015, the quasars were 
taken at the distances z < 0.3, and BL Lacs were se- 
lected with luminosities m < 18 as in 7 . According to 
Fig. [3[ there is no significant deviation of the number 
of coincidences from isotropic expectation in the case of 
Lacertae and quasars. Variation in the redshift ranges 
for quasars does not reveal any significant deviation. 
Any significant correlation is also absent for HP objects 
and BL+HP objects from [8]. 

Correlations with AGNs are found in the PAO and 
Yakutsk array data at almost the same angular dis- 
tances (3.1° and 3°, respectively). It is unclear whether 
this is a random coincidence. The error in arrival an- 
gles of the showers with axes inside and outside the ar- 
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ray area is estimated in p3] as 3° and 5°, respectively. 
At the same time, as the energy increases, the angular 
error should decrease because the number of array sta- 
tions fired in a shower increases with energy. Therefore, 
the actual accuracy of UHECR angles for the Yakutsk 
array can be equal to 3° or better. 

Approximate coincidence in the threshold energies 
at which the maximum correlation is observed by the 
two arrays may be accidental. A comparison shows the 
similarity in the shapes of the energy spectra measured 
by giant EAS arrays: if the correction factors are applied 
to the energy estimates of the primary particles induc- 
ing EASs, the spectra almost coincide both in shape 
and in intensity [20] . Systematic energy corrections in- 
troduced in this way for the Pierre Auger observatory 
and Yakutsk EAS array can differ from each other by a 
factor of 1.5 — 2. Hence, the threshold energy at which 
the maximum correlation is reached for the PAO data 
would be from 90 to 120 EeV if the assumption about 
correction factors is valid. 

Another part of the PAC hypothesis is that cosmic 
rays correlate with the directions to the AGNs at dis- 
tances z < 0.018 due to the Greisen-Zatsepin-Kuzmin 
effect 21;, which strongly suppress the flux of cosmic 
rays with energies E > 60 EeV from cosmological dis- 
tances. In order to verify this effect with the Yakutsk 
array data, the ratio of the observed number of coinci- 
dences of cosmic rays (E > 60 EeV) with AGNs in vari- 
ous redshift bins to the number of random coincidences 
expected for the isotropic distribution is used (Fig. 2]). 
Indeed, a significant correlation of UHECRs with AGNs 
is found only in the interval z £ (0.001,0.015). In all 
other redshift bins, the observed number of coincidences 
is equal (within errors) to the number expected in the 
isotropic case. Hence, this can be considered as one 
of the independent evidences of the Greisen-Zatsepin- 
Kuzmin effect. 

It should be noted that the Yakutsk data, as in the 
case of the PAO data, possibly correlate not only with 
AGNs, but also with other astrophysical objects with 
a similar spatial distribution that are not considered in 
this work. Moreover, Gorbunov et al. [22] stated con- 
trary to the PAC hypothesis: the conclusion that the 
most part of UHECRs are protons originating in nearby 
extragalactic sources (AGNs) can be rejected at 99% 
CL. Instead, they explained the data P2] by the exis- 
tence of a bright source in the direction of the Centaurus 
Supercluster. Another interpretation of these data was 
proposed by Wibig and Wolfcndale [23] , who stated that 
cosmic rays are nuclei with < InA > = 2.2. ±0.8 that are 
generated in nearby (about tens of Mpc) radio galaxies. 
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Fig[5] Ratio of UHECR coincidences with different 
classes of AGNs to the expected number. 



Analysis of the muon detectors data of the Yakutsk 
array [24] shows that the data at energies above 10 
EeV can be explained within the framework of the two- 
component model of the cosmic ray composition consist- 
ing of protons and a considerable fraction of heavy nu- 
clei. This is hardly consistent with the correlation of the 
arrival directions of such particles within 3° with AGNs. 
'Clusters' of EASs associated with clusters of (radio) 
galaxies are hardly seen in the data of the Yakutsk ar- 
ray in contrast to the PAO data. 

If the hypothesis that active galactic nuclei are 
sources of UHECRs is accepted, the following ques- 
tion arises: Whether all classes of the objects belong- 
ing to active galaxies are sources or, for example, only 
Seifert galaxies generate cosmic rays, as was assumed 
by Uryson [llj . To answer this question, one can an- 
alyze a correlation of UHECR arrival directions with 
each class of the objects separately. Figure [5] shows the 
results of such an analysis of the data. Here, AGNs 
are divided into four classes according to the proposal 
of VeronCetty and Veron [8]: i) SI, Seifert galaxies of 
the first type with broad Balmer lines; ii) S2, Seifert 
galaxies of the second type; iii) S3, so-called LINERs 
(low-ionization nuclear emission- line regions) , which are 
galaxies with weak nuclear emission lines; and iv) H2, 
galaxies whose spectrum of nuclear emission lines is sim- 
ilar to that of nebulae ionized by hot stars. The redshift 
boundaries giving the maximum correlation with the 
Yakutsk data are selected for each class of the objects: 
z < 0.015 for SI and S3, z < 0.016 for S2, and z < 0.024 
for H2. 

Correlations of the Seifert galaxies of the first type 
and LINERs with cosmic rays do not exceed correlation 
expected for the isotropic distribution. Only S2 and H2 
objects correlate with UHECRs (with the maxima at 
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E > 50 and E > 70 EeV, respectively). The excess of 
the observed number of coincidences over the number 
expected for the isotropic distribution is 4.7(7 and 3.7er 
in the maxima for S2 and H2, respectively. Therefore, 
possible sources of UHECRs are Seifert galaxies of the 
second type and/or H2 objects at distances less than 
100 Mpc. 

This work is supported by the Russian Foundation 
for Basic Research (project no. 06-02-16973) and by 
the Council of the President of the Russian Federation 
for Support of Young Scientists and Leading Scientific 
Schools (project no. NSh-7514.2006.2). 

Translated by R. Tyapaev. 
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